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Abstract 

The functional determinant of an elliptic operator with positive, discrete spec- 
trum may be defined as where Z{s), the zeta function, is the sum J2'^ 
analytically continued to s around the origin. In this paper Z'{0) is calculated for 
the Laplace operator with Dirichlet boundary conditions inside polygons and sim- 
plices with the topology of a disc in the Euclidean plane. The domains we consider 
are hence piece-wise flat with corners on the boundary and in the interior. Our 
results are complementary to earlier investigations of the determinants on smooth 
surfaces with smooth boundaries. The computation takes the form of a variation of 
the shape of the domains, which is chosen such that the coordinate transformations 
are conformal everywhere but at the corners. The contributions to the variation of 
Z'(0) in simplices then come almost exclusively from terms with singular support at 
the corners: in polygons there are no contributions but from the corners. We have 
explicit closed integrated expressions for triangles and regular polygons. Among 
these, there are five special cases (three triangles, the square and the circular disc), 
where the Z'{0) are known by other means. One special case fixes an integration 
constant, and the others provide four independent analytical checks. 
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1 Introduction 



One of the basic integrals that arises in many parts in physics is 

n 7^^"^^^^ = (DetA)-^^ (1) 

where A is a real, symmetric matrix with posive eigenvalues. For instance, let (|l]) describe 
the integration of fluctuations around a classical solution in quantum mechanics, where 
the Lagrangian has been expanded up to second order. The dimension of A is then infinite, 
and (|l]) is divergent on both sides. The equation is therefore undefined as it stands. As 
a basic example take a one-dimensional harmonic potential, and the fluctuations in the 
time interval < r < L. Then 

xAx = / dT{drxf + u; V, (2) 
Jo 

and we have Dirichlet boundary conditions for x at r = and t = L. The most straight- 
forward way, is then to go back to the Gaussian integral (|l|), reintroduce a cut-off e in 
space, and modify the integration measure depending on the cut-off so that the limit 
when e goes to zero is flnite. In quantum mechanics this is feasible: changing (27r)~2 to 
(27re)~2, and including one more factor (27re)~2, turns (|T]) to a discrete approximation to 
Feynman's sum over paths, which in the limit gives 

, . -.r^ .nx 1 ,27r sinhftuL) , i . , 

(Regularized [Det A] )-5 = ( i)-2, (3) 

and this is the correct expression in the Greens function. 

We may also observe that the eigenvalues of A are ^^j^ + u"^, and the determinant is 
then formally 

detA~n(^ + ^') (4) 

n=l ^ 

One way to regularize the determinant is to introduce a cut-off A in the product @, 
check that in the limit of large A the result separates into one flnite factor and one factor 
divergent with A, and keep the flnite factor as the regularized result. For (^, this gives 
the same result as (PD[P^. 

A regularization can also be found from zeta function of the operator; 

Za{s)=Y.K' (5) 

which converges when the real part of s is large enough. When this function can be 
analytically continued to be regular in a neighbourhood of the origin, then 

Regularized [Det A] = e~^^(o) (6) 
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For (^), this again gives the same result as (^. The zeta function method was first 
introduced in the context of regularizing expressions like (|I]) by Hawking|^, to study 
fluctuating fields in a background of curved space. 

In dimension higher than one, it is not evident that all regularizations of the deter- 
minant give the same result. In this paper we will compute Z'j^{0), with A the Laplace 
operator with Dirichlet boundary conditions in two-dimensional simplicial domains. In 
the language of field theory this is the gaussian model, a free massless theory, albeit with 
unusual and obstructing boundary conditions. Even so, the direct limit is far from trivial, 
and the most relevant result we are aware of, works only for lattice Laplacians, discretized 
on ectangular (M x A^) domainsj^: 

detA^ 2te^(l + ^^2)-^ (MiV)-ir^(g)(^)i (7) 

Here, in lattice units, MN is the area, 2{M+N) is the length of the boundary, q = e~^'^^ is 
the modular parameter, G is Catalan's constant and 7]{q) is the modular form of Dedekind. 
There are now no less than three terms separately diverging with the size of the lattice. 
If, with hindsight, we use that for rectangular domains Za{0) = |, we can rewrite (|^) in 
terms of an explicit lattice spacing a: 

Area Length 

det A ~ /i7^/i7^(a=')^-^(o)e-^-^(°)'°s^"'^"-^ (8) 
where are the various remaining terms in which agreeQ with e~^^*^°''[0](see Ap- 



pendix 11). If nothing else, it seems likely that a discretization on a rectangular grid, of 
a domain which is not itself of rectangular shape, will give rise to oscillating terms in the 
cut-off. If d^) is to be generally valid in two dimensions, it can probably only be of a 
smoothened discretized determinant, where the smoothening goes over cut-off scales. As- 
suming that this can be done, and considering that the area, the length of the boundary, 
and Za{0) are all integrals of local distributions (see section 2), it is possible to introduce 
local cut-off dependent counter-terms, such that the the finite remaining piece is e~^'A^^\ 

It therefore at least makes sense to define the regularized determinant to be 6"^^^°-*, 
and this is the view we take in the rest of this paper. We will use the notation Z'j^{0) for 
our generic case: the zeta function of the Laplacian with Dirichlet boundary conditions in 
a simplicial domain D, with the topology of a disc. We will freely change the index of the 
zeta function to denote various special cases, and even contributions to the regularized 
determinant from some parts of the domain. 

If we look for physical relevance, we must have fiuctuating geometry, as otherwise all 
would dissappear in an overall normalization. It is a quite old idea that the elementary 
excitations of non-Abelian gauge theories are string-like object s[pT|, \n\ , |T|. In lattice 
gauge theories, the statistical weight of a Wilson loop, when a quark and an anti-quark 
are taken apart for some time, is the area of the smallest area delimited by the loop. A 
phenomenological model of the excitations was proposed [|1^], where the statistical weight 



^up to a constant factor 2^4 
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of a surface is its area, and the path integral goes over imbeddings in d-dimensional 
external space ("x'^"), and over internal two-dimensional geometry ("(7"^"): 

Our computations are relevant to a part of the investigations of (||). With g""^ fixed, 
the integration over is just a quadratic integral with Dirichlet boundary conditions 
like ([l|). Our calculations hence gives the finite piece of this determinant. As is well 
known, reparametrization invariance of the action in (|^), gives rise to Faddeev-Popov 
determinants, which for smooth surfaces turn out to be determinants of Laplacians acting 
on vector fields, with modified Dirichlet boundary conditions. We have not investigated 
these determinants, and we are not quite sure they are relevant when we consider piece- 
wise flat surfaces with sharp corners. At least in a class of simplices with fixed number of 
corners, such a surface is its own model, and we would not have any more reparametriza- 
tion invariance. In this respect, an approach closer to the simplicial discretization of (P) 
would seem to be more appropriate 0]. 

Let us now see why it may be interesting to investigate the determinants on simpli- 
cial disc-like domains. Smooth disc-like manifolds with smooth boundary can always be 
mapped conformally onto one another. If we denote the conformal factor a{x), the base 
metric and curvature by g and R, and the base geodetic curvature of the boundary by k, 
a celebrated result [l[] says that 

Z^,(0) = Z'^{0) -—I dsh-da+ — f dska 

4:71 JdD 67r JdD 

I d^zJ~g[r''daad^a + Ra] (10) 
lZtx Jd * 

The integration constant can be computed from the upper half sphere[^]. If we disregard 
the boundary terms, the equations of motion of are Liouville's equation, and we refer 
to the combined action in (|TD|) as the Liouville action. 

In ordinary Feynman path integrals (|I|) the typical path is quite rough, i.e. a nowhere 
smooth random walk. It seems likely that the typical surface that enters in ^ is also 
quite rough. On two-dimensional smooth surfaces, one can open up a corner with angle 
27ra (or vra at the boundary) with a coordinate transformation which is conformal and 
regular everywhere but at the corner, where it instead has a logarithmic singularity. The 
kinetic energy term in (0) will then be logarithmically divergent at the corner. In other 
words, the action for smooth (and conformal) coordinate changes is infinite for these 
transformations. One possible procedure is then to introduce a cut-off vq at a corner, 
and simply remove the quantities diverging when tq goes to zero|^. This is obviously 
dangerous; if one is not careful one easily gets a spurious finite piece from a logarithmic 
divergence, and Z'j^{0) is a well-defined mathematical object, which has a value, the surface 
being smooth or not. The correct interpretation is that a conformal transformation does 
not change Zoi^O), but a non-conformal one does. In fact, for piece-wise fiat surfaces, 
Zd^O) can be written as a sum over a rational functions of the opening angles of the 
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corners [see section 2, Appendix 6). A typical non-conformal transformation will change 
the opening angles by a shear, so it will change Z£)(0). The infinity of (|I^) under non- 
conformal transformations is therefore a mirror of actually Z£){0) changing. 

Our computation give some more explicit results on determinants, to which one does 
not have access from smooth models. This can have some mathematical interest by itself. 
More speculatively, it is possible that a definition of determinants by a precise calculation 
of Z'jj{0), may yield a better regularization of (|]) than does (p!0| ) and its Faddeev-Popov 
ghosts. Certainly, such a result would go far beyond what is actually done here: we have 
not begun to address a computation of a sum over surfaces as in (|^). 

The conclusions of this paper can now be stated as the following propositions: 

Proposition 1. Domains with disc-like topology may be mapped to the upper half com- 
plex plane, the boundary being mapped to the real axis, and the corners being mapped to 
branch-points ujj . If z is the coordinate in the domain, the transformation satisfies 

^ = 0i.^(^-^.^"^ (11) 

and the representation is determined by the angles and the lengths of the sides, up to a 
rational fractional transformation of the upper half plane. Choosing one parametrization 
the normalized area of the simplex is: 

r duj A dui T-r , ,9^ 9 / N 

and the variation of Z']j{0) under a general shear and dilatation can be written; 

S[ZDmArea=A] = ^[^0(0)] + ^[Zd{0) logA] - 6[Zn{0) log Area^] (13) 
where then value of the zeta function at the origin is; 

interior corners ^ boundary corners 

and the opening angles are written 2Ttaj in the interior, and naj on the boundary. 

Proposition 2. When the area is chosen Areao, the difference of ^^(0) between two 
simplices, that differ by an infinitessimal transformation, which is regular and confor- 
mal everywhere except at the corners, can be written as a sum over quantities with point 
mass support at the corners, and a line density on the boundary. The integral over of the 
boundary term is identically zero for polygons, but gives 

-47r E ^^i^ (15) 

interior corners 
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if there are corners in the interior. 

Proposition 3. The contribution to 5Z^(0) from the point mass in one corner c, is 
a sum 

(5Z^(0)|c = 5Z'^^{Q) + Z«^(0) • 5[other corners] (16) 

where (0) is exclusively determined locally at the corner, -^ac(O) 'the contribution to 
Zd(0) from the corner, and the influence from the other corners depend on the opening 
angles of these, and on the lengths of the sides: 

S(oj LO ) 

Zad^) [other corners] = Zq^(0)[^ 5q;c' log |a;c' —^d^ — (1 — q;c')[ ' ^ + c.c]] (17) 



Proposition 4. The strictly local contributions may be given in integrated form, and 
are for a corner on the boundary: 

Z'aiO) = :j^(- - «)(7 - log2) - 1(1 + 3 + a) loga + J(a), (18) 
12 a 12 a 

and for a corner in the interior; 

Z'M = - a)(7 - log2) - 1(1 + a) loga + 2J{a), (19) 

where the term J has the integral representation; 

J{a) = r ^[;l(coth(f ) - acoth(|)) - -1(1 - a)]dx. (20) 
Jo — 1 2x 2a 2 L2 a 

We have here included an integration constant in our definition of J. 

Special situations are quite important to us, as they provide necessary checks. We there- 
fore hst the them separately: 

Proposition 5. For triangles one may choose a parametrization where the branchpoints 
lie fixed in 0, 1 and oo. The normal area of a triangle with angles nai, na2, na^ is then 

A ( ^ ^ r(ai)r(a2)r(a3) /^^v 

Area{ai, as, as = 77 vFTl vFri V (^1) 

2 r(l - Q;i)r(l - a2)r(l - a^) 

and all terms but the strictly local in the variation o/Z^(0) vanish. The determinant of 
a triangular domain is thus 

^t(0)= E ^;.(0) (22) 

p=l,2,3 

with Z'^{0) as in Proposition 4- 
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Proposition 6. For regular polygons with n corners one can choose a parametrization by 
mapping to the unit disc, where the corners are regularly spaced on the unit circle. The 
radius of the circumscribed circle is with this parametrization: 

r(i + i)r(i-^) , , 

= rn -1) ^^^^ 

and the determinant in a regular polygon with radius of circumscribed circle R is; 

Z'p„ (0) = Zp^ (0) log ^ - — ^ log n + nZ[_.{0) (24) 
where Zp^{0) = q^-^- Specializing to n = 2 we obtain the determinant in a square: 

Z'pM = \ log Area + ^ log ^ + ^ log tt + ^ log 2 (25) 
Taking the limit as n ^ oo we obtain the determinant in a disc: 

Z'pjO) = i logi? + ^ + 1 logvr + i log2 + 2C'(-1) (26) 



The paper is organized as follows: standard results on heat kernels and zeta functions 
are summarized in section 2. In section 3 a representation is found for 6Z'{0) which 
treats the corners and the rest in different ways; this is the main idea. In both cases 
the representation is in terms of the short-times asymptotics of the heat kernel. In the 
interior this representation is standard, and gives the variational form of ([T0|). That the 
underlying space is fiat, and with disc-like topology, actually simplifies things so much 
that this contribution gives only the simple sum in Proposition 1. In the corners the 
appropriate asymptotics is the Sommerfeldt heat kernel in an infinite sector. Most of the 
derivations are by calculation, and most are elementary (although sometimes cumbersome) 
We have chosen to present those in appendices, largely in a self-contained way, but without 
any reference to physical arguments. The appendices can therefore be read more or less 
independently from the rest of the paper. Section 4 summarizes the results. 



2 Heat kernel and Zeta function phenomenology 

The diagonal elements of the heat kernel on smooth manifolds admit an asymptotic ex- 
pansion for short times. In two dimensions the expansion goes as 

Kj,{x, X, t) ~ + -i^ + Co{x) + ... (27) 
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McKean and Singer |T5[ proved the existence of this expansion for smooth manifolds up to 



terms t, and for manifolds with smooth boundary up to t^. In the interior the coefficients 
Cj are polynomials in the curvature tensor, and on the boundary they are line distributions, 
with weight depending on the curvature of the boundary. For polygonal domains in the 
plane, Kacfl^ proved the expansion up to co(x), which in this case is made up of point 



masses at the corners. Since the curvature of a polygon can be said to be concentrated 
at the corners, the result is in some sense natural, and Cq{x) for a smooth boundary can 
indeed be found as the limiting case of polygons with angles coming closer and closer to tt. 
However, the converse is not true, that is, if one wants to compute Z£,{^) for a manifold 
with tips and corners, it is not possible to use a smoothened approximation. 
The integrated form of ( P7D is 



/ dx^Koix, X, t) = &D{t) ~ ^ + ^ + Co + . . . (2^ 

JD t +2 



The zeta function is defined by a Mellin transform of the trace of the heat kernel as 



1 

Zn{s) = — - / dW-'eoit) (29) 
i [s) Jo 

The integral only converges in the lower limit if > — ctj for all terms Cc-t"' in the 
heat kernel, i.e. 3fJs > 1 in two dimensions. One can get around the pole by a partial 
integration, where the boundary terms vanish for 3fJs > 1: 

/ dtf-^eoit) = 

Jo 



r{s) 

■[( Tt^Dimo-T^J dt ~[tQD{t)] (30) 

s — 1 i[s) Jo s — lot 



Vis 



The second integral now gives an analytic continuation down to > |. Further on the 
zeta function has poles at s = — ««, with residue -^rH, except at zero and the negative 
integers, where the inverse gamma function has zeros, which gives finite values Z{—n) = 
(-)V_„n![H. 



For positive n greater than 1, a method that goes back to Euler[l£] gives a represen- 



tation of the zeta function as a convolution of electrostatic Greens functions: 

Znin + 1) = rr^(^)"+^ = J dx J dyi . . . J dynGoix, yi)GD{yi, Z/2) • • • Goiyn, x) (31) 

This method can be extended to calculate the finite part of Zd{s) at s = 10. Hence, if 
the expansion (p8[) goes on indefinitely and one knows the Greens function, one may in 
principle calculate the zeta function at the integers, and the residues at all the poles. But if 
this is true, and the zeta function does not grow too fast at infinity, then the zeta function 
is completely determined, which (somewhat indirectly) determines the eigenvalues. Hence 
all the information about the eigenvalues, and on all quantities depending on them, is 
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then contained in the Green's function and the asymptotic expansion of the heat kernel 
for short times. 

It is natural to consider what we will call a zeta function density: 

1 r°° 

Zd{x,x,s) = -— dtf-'KD{x,x,t) (32) 

L{S) Jo 

For < 1 we define the zeta function density by analytic continuation from (^21). At 
zero, the zeta function density is determined by the asymptotics of the heat kernel for 
short times as in (pTl) : 

Zd{x,x, 0) = Co(x) (33) 

We may just as well consider the density associated to the derivative of the zeta function 
at the origin, but this quantity will depend also on the heat kernel for long times. These 
questions are addressed in more detail in Appendix 1. 

One could also look at off-diagonal elements in (|32D , or equivalently, the operator 
(z^)^- Since the off-diagonal Mellin transforms are finite for all s, the analytical contin- 
uation brings no off-diagonal subtractions, and effectively the diagonal and off-diagonal 
elements are treated differently for 5Rs < 1. It would be interesting to have a proper 
regularization of off-diagonal terms, so that one might consider Z'^{0) as the trace of 
a (regularized) operator Q ~ — log(— A), and in terms of which a regularized Laplace 
operator would be e~'^ (in the operator sense). In this paper we will only consider the 
diagonal elements, and we therefore refrain from calling the zeta function density the 
matrix elements of an operator when 3?s < 1. 

Let us consider a dilation x — Ax of the domain. That changes the eigenvalues 
in a simple way: En — > X~'^En, and the zeta function accordingly also changes simply: 
Zxd{s) A^*Zd(s). For the quantities arond the origin this implies: 

Zxd{0) = Zd{0) Z'^^iO) = Zn'iO) + \ogX^ZniO) (34) 

The invariance of Z£,{0) under dilations is a special case of invariance under regular 
conformal transformations. The change in Z'^{0) logarithmically proportional to the 
variation of the area, is contained in Proposition 1, and was already noticeable in the 
asymptotic expansion of the discretized determinant in rectangular domains @. 

3 A variational formula 

In this section we will write down a variational formula. It can be formulated either 
in terms of the asymptotics of the heat kernel for short times, or in terms of the zeta 
function density from section 2. Let us first motivate why we give two computations, that 
essentially only differ in that they work on opposite sides of the Mellin transform: the 
approach using the heat kernel is more traditional, but requires more delicate treatment 
to separate out the finite piece. Note that our formulae also differ by (a fairly simple 
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term) from a definition of the regularized determinant often used in field tlieory[0, [T] 
The reason for this discrepancy is that we also consider transformations that are not 
conformal, and here the additional term matters. The approach using the zeta function is 
computationally somewhat simpler, once the notion of a zeta function density is admitted. 

Let now D and D' = D + 6D be two domains that differ infinitesimally. Then, for 
sRs > 1, 



oo 



Zd'{s) - Zn{s) = - E-'] = E(^)^n^ (35) 

and this definition is analytically continued to lower values of s. We can consider {-E^} to 
be the eigenvalues of a modified operator (—A — 5A) in the domain D, so that we write 

Zd.{s) - Zj,{s) = sTr,[{5I\)ZD{x,x,s + 1)] (36) 

The variations we will consider map one simplex on another. They therefore leave the 
following sum invariant: 

S{ai) = + E + _ _ E (2 - 2a,) = 2 (37) 

boundary corners interior corners 

where the opening angle of a corner on the boundary (in the interior) is vra {2T^a). The 
analytical expression for our variation will therefore only be determined up to the variation 
of a (s-dependent) function of 5(aj), that is, up to terms linear in the variations of the 
angles. 

If we want to evaluate the variation of the derivative at zero we can do it as 

^5Zo{s)\s=o = \im-Tr^[{ZD'ix,x,s) - ZDix,x,s)) - {Zn'{x,x,0) - Zoix,x,0))] (38) 
as s^o s 

which can be rewritten as 

-^6ZDis)\s=o = Fimtes^o^Tr^[ZD'ix, X, s) - ZDix,x,s)] (39) 



In general there will be a pole at the origin in s in (|39D, with a residue equal to the 
variation of Zd{0). Combining (^) and (|36|) we have 

-^6Znis)\s=o = Fimtes^oTr^[{5A)Znix,x,s + l)] (40) 

It is now convenient to choose the variation in a particular way. Think first of the Laplace- 
Beltrami operator in curved space, with metric tensor g"-^: A = g~2dag°'''g^db. In two 
dimensions we can choose the particular form gab = e'^'^^^^gab, and for a smooth surface 
with disc topology we can take g"''' to be the standard flat metric. In this coordinate 
system, the Laplacian has the following form A = e~^'^(9^^. A conformal variation of the 
Laplace-Bertrami operator is then 6 A = {—26a) A, and the final variational formula is 

-^SZois)\s=o = Fimtes^oTr^[25a{x)ZDix,x,s)] (41) 
as 
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The formula (^) holds also when we include corners with varying opening angles, the 
only new effect then being that prefactor weight Sa has a logarithmic singularity at the 
corner. 

We will show in Appendix 3 that in fiat domains far from the boundary, the density 
Zd{x, X, s) vanishes at the origin. Hence the interior of fiat domains, far from the corners, 
give no contribution at all to Z'^{0). This is in agreement with (P), since our base space is 
fiat (hence R is zero), and the kinetic energy term of a regular conformal variation can only 
give a boundary term. At straight boundaries, there is an undetermined contribution to 
Z^(0), proportional to the boundary length. The contribution to the variation, 6Z'j^{0), 
is however fully determined following (|4l|) , and picks out the normal derivative of the 
variation 6cr{x). Here again we have exact agreement with (p!0|). 



For simplices, the boundary term integrates trivially, and we are left with only the 
corner contributions as the important parts. Considering that the zeta function density 
Zn{x,x,0), has components with point mass support at the corners, it may be surmised 
that integrated against the logarithm of the distance from the tip, it gives an infinite 
contribution to (|^). This is indeed the case, that infinity is precisely the variation of 
Zd{0), which we subtract by taking the finite part as s goes to zero. It is now important 
to realise that we can choose small areas around the corners, with some radius tq, which 
we can let tend to zero at the end of the calculation. It therefore does not matter if we 
compute the answer up to terms proportional to tq or Tq, since these will eventually drop 
out. This means that we can actually substitute the true (unknown) Z£,{x,x, s), with a 
sufficiently good approximation, computed from the Mellin transform of an approximation 
to the true heat kernel valid for short times. The technique for doing this is explained 
in Appendix 5, and the approximation to the heat kernel - the Sommerfeldt heat kernel 
in a sector - is described in Appendix 4. The calculations of the contributions to Z'^{0) 
from the corner is then done in Appendix 6. 

Let us now derive an alternative formula to (^TJ) using only the heat kernel for short 
times. We begin with the following equality (derived in Appendix 1): 

poo (If 

Z'^{0) = -fTr^Zoix, X, 0) + Finite.^oTr, / -/^^(x, x, t) (42) 

J t t 

The heat kernel can be expanded in a complete set of states: 

KD{x,x,t) = J2\M^)\^e-^^' (43) 

V 

and the variation of (^2|) will be 

roc 

SZ'o{0) = -f6Tr^ZD{x,x,0) - Finitee^o / dtY, < v\{26aA)\v > e"^"* 

= •ySTr^ZD^x, x, 0) + Fimte^^oTrx[26a{x)K£){x, x, e)] (44) 

When we are in the interior, or at a straight boundary, the zeta function density ZdIx, x, 0) 
is zero. The remaining trace over the heat kernel in (HI) then agrees with a well-known 
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definition of tlie regularized determinant in field theory, and the short-time expansion of 
the heat kernel can be done by standard means p!7|, We then find the variational form 
of ([T0|). In piece-wise fiat simplices there are further simplifications, such that the only 
non-zero term is the integral over the boundary of the normal derivative of 5a in (|TUp. 

Close to the corners, we cannot however use the same short-time expansion of the 
heat kernel. Again we have to turn to the Sommerfeldt heat kernel, and carefully extract 
the constant piece of the trace of the heat kernel for short times. Finally, we must take 
into consideration that the variation of Zi?(0) is not zero at the corner, and add the first 
term on the right-hand side of The computation following these lines is done in 

Appendix 7. 



4 Calculations & Conclusions 

In this section we describe the computations in the appendices, and discuss the results. 

Appendix 1 is an expansion on Section 2. We establish that the zeta function density 
of the derivative at zero is well-defined, but depends on the heat kernel for long times. 
An alternative definition of the variation of the regularized determinant in terms on the 
heat kernel for short times is derived. Appendix 2 presents a general parametrization of 
variation of simplices, and computes the contribution from corners in the interior as in 
Proposition 2. In Appendix 3 we treat in some detail the case of a straight boundary. 
We derive that it cannot give a a term in Z'j-,{x,x,0) more singular than a line density, 
and hence a contribution to Z'j^{0) proportional to boundary length. We also derive the 
contribution from the line integral on the boundary to the variation of Z'^{0). Appendix 4 
presents the Sommerfeldt heat kernel in a sector, and Appendix 5 states an elementary 
integral, that is convenient when one wants to compute the Mellin transform using the 
Sommerfeldt heat kernel. 

In Appendix 6 we compute the contribution to Ze,{0), and the strictly local contribu- 
tion to Z'j^{0) from a corner with opening angle a, using the approach of the zeta function 
density. We call these quantities Za{0) and Z'^{0). We perform the computations for both 
corners on the boundary and corners in the interior. In Appendix 7 we do the same thing, 
using the slightly more involved approach of the analysis of the short-time properties of 
the Sommerfeldt heat kernel. In Appendix 8 we investigate Z'^{0) in the special cases, 
where the opening angle is -. In Appendix 9 we investigate the asymptotic behaviour of 
Z'^{0) as a is large, small or close to one. 

In Appendix 10 we treat a few special cases using our general variational formula. The 
additional computations done here essentially boil down to taking care of the expression in 
Proposition 3, that has support at one corner, but depends on the angles of all the corners, 
and on the lengths of the sides. In the very special case of triangles, this additional term 
is absent, but in polygons it is important. In Appendix 11 we state summarily the cases 
we know to compare with, where the Z'j^{0) are known, either because the eigenvalues 
are known, or because these domains are conformal images of some domain where the 
eigenvalues are known. This is important information: it fixes an integration constant 
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that we cannot reach; and it provides much needed analytical checks. 

In conclusion we have shown that the determinants on piece-wise smooth surfaces are 
different from the determinants on smooth surfaces. The formula that describe variations 
of determinants on smooth surfaces (^) give infinite results for variations that change 
the opening angles of the corners. This infinity is not real, it is only an imprint of 
actually Zij{0) changing. We have parametrized the variations of the domains such they 
are conformal everywhere outside the corners, and separated the contributions to the 
variation of Z'^{0) into four parts: one (Proposition 1), that expresses the variation of 
Z£){0) times the logarithm of the area; one (Proposition 2), which is an integral over the 
boundary of the normal derivative of the conformal parameter; one (Proposition 3) which 
has local support at the corners, but depend on the opening angles of all the corners and 
on the lengths of the sides; and one (Proposition 4), which is determined strictly locally 
at the corners. 

The first part amounts to determining how the normalized area changes with the 
angles and the sides. We can trivially integrate the second part, which is identical to 
one of the terms in the variation of determinants on smooth surfaces (p!0|). We have 
integrated the fourth part, both for corners on the boundary and in the interior. For 
general opening angles this part has an integral representation ( |A6.29| , |A6.30D , which for 
the special opening angles a = ^, can be resolved into a finite sum ( A8.11|) . 

We have been able to parametrize the area and integrate the third part, for the special 
cases of triangles (where the third part vanishes) and a class of polygons, which includes 
certain interpolations between pairs of regular polygons (Appendix 10). In general, we 
have not been able to integrate the third part, which therefore has to be left in variational 
form (Proposition 3). 

The wider applicability of our approach evidently depends on better handling of the 
parts we have not integrated. We do not expect that one will in general be able even 
to express the normalized area in closed form, as a function of the angles and the sides. 
Perhaps one may however find a limiting form, which is valid for small opening angles. 
Considered as an action for simplicial surfaces, it is possible that the full expression of 
Z^(0) strongly damps out corners with small opening angles. As is well known, most 
investigations of sums over random surfaces in physical dimensions lead to very rough 
surfaces, dominated by long sharp spikes. If the full expression of Z'j-){0) as a function of 
one opening angle has a singularity at the origin - with the right sign of the prefactor - 
then surfaces with spikes will be more strongly damped by the action Z^(0), then by any 
finite expansion in local curvature. We hope to return to questions in this direction in 
the future. 



5 Acknowledgements 

E.A. wishes to thank Claude Itzykson for initially encouraging him to study these prob- 
lems, and Service de Physique Theorique at Saclay for hospitality in 1986-87. We thank 
Jean-Marc Luck for communicating numerical results on the regularized determinants in 



13 



triangles, that proved a valuable check. This work was supported by the Swedish Natural 
Science Research Council under contracts U-FR-1778-101, S-FO 1778-302 and F-FU 
8230-306, and by the Goran Gustavsson Foundation. 



14 



Al. Local and nonlocal quantities 

In this appendix we collect some results on the zeta function density and compare it to 
other approaches in the literature. An important aspect is locality: a Mellin transform 
of a heat kernel, at a point, integrated over some subdomain, or over the whole domain, 
is said to be local, if it is determined only by the expansion of the heat kernel for short 
times. 

We start with the definition 



1 p'^ 

Zd{x,x,s) = -— dW-^KD{x,x,t) (Al.l) 

1 S JO 



It is useful to explicitly analytically continue the integral beyond the poles at s = 1, s = | 
and also beyond s = 0. We will then get out a pole in s at s = from the integral, which 
will combine with the gamma-function in front to give a regular function. We thus have 
a representation of the zeta function density around s = in terms of a regular prefactor 
and a convergent integral: 



1 r°° d id id 

^°'"-^--^)= (.-l)(.-lKw /o ^f-\t^giPgitKn(^,-M (A1.2) 



1 r"^ ... d ,.1 d ,.1 d 

Using the asymptotic expansions as t is close to zero: 



C\(x) Ci[x) 1 
KD{x,x,t) ~ + h co{x) + C_i{x)t2 + ... 

d ici{x) 3 1 

—tKD{x,x,t) ~ h Co(x) + -C_i(x)t2 + . . . 

1 1 d 1 3 1 

tH-Q^[t^Q^tKD{x,x,t)] ~ -co(x) + -c„i (a;)t^ + . . . 

we may evaluate as follows 

ResZ{x,x, s)\s=i = ci(x) 

ResZ{x,x, s)\g^i = 7r^ci(x) 
Z{x, X, 0) = co(x) 

We may also evaluate the derivative with respect to s at the origin: 

/■oo Q ^ Q ^ Q 

Z^(x,x,0) = (3 + 7)co(x) + 2y^ dt\ogt—p — p—tKD{x,xM (Al.3) 

It is quite clear that Z'^{x,x,0) is a quantity which depends also on the heat kernel at 
large times. 

An alternative formula is obtained by taking the integral from e to infinity, where e 
is some small positive number. As the integrand in [A1.3| behaves as t~2 logt for small t, 
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the integral is convergent, and the hmit as e goes to zero is harmless. Partial integrations 
will however bring in divergent terms from the lower boundary. A short calculation gives 



Z'jj{x, X, 0) = (3 + ^)co{x) + 

1 f°° dt 

(co(x) loge — 2ci (x)e^2 — ci(x)e^"'^ — 3co(x)) + / — KD{x,x,t) 



which can be simplified to 

poo (^f 

Z'^{x, X, 0) = 7Co(x) + Finitee^o / —Kd{x, x, t) 



(A1.4) 



(A1.5) 



Let Kd be an asymptotic approximation to K^,, valid for short times. We may then write 
the integral in |A1.5| as 



1 dt 



t 



1 dt 



Finitee^o / — KD{x,x,t) + / — [KD{x,x,t) — KD{x,x,t)] + / — KD{x,x,t){Al.6) 



t 



dt 



1 t 



The second and third terms are finite, if the asymptotic approximation is sufficiently 
good. Hence they can only contribute to Z'j^{x,x,0) as a smooth density. The only term 
which can give a distribution contributions, the space integral of which contains a finite 
piece even if taken over a vanishingly small area, is the first in [A1.6| , and the term 7Co(x). 
Both only depend on the heat kernel for short times. 

When we integrate [A1.5| over the whole domain, we obtain 



/CO df 
jKDix,X,t) 



(A1.7) 



A2. Variation and parametrization 

We normalise the areas by choosing the simplest form of the Schwarz-Cristoffel transfor- 
mations that map the disc-like domains to the upper half complex plane. Let z be the 
variable in the domain, and u the variable in the upper half plane. Then we have: 

^ = ^^Ylico - co.y--^ (A2.1) 
j 

The vertices of the domain map to the branch points ujj, which are determined up to 
a rational fractional transformation of the upper half plane. This redundancy may be 
eliminated (up to permutations) by fixing the positions of the images of three vertices. 
An overall phase factor 0d is determined by the orientation of the domain D in the z- 
plane. For the simplest case, triangles, the three branch points may be put in 0, 1 and 
oo in the t^-plane, and the side of the triangle between the vertices that are mapped to 
and 1 can be placed parallel to the real axis in the 2;-plane. We then have the familiar 
formula: 

= cji-"o(i (A2.2) 

dz 
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A variation of the shape of the domain, brings a simultaneous variation of all the 
angles, and of the positions of the branch points: the expression for the variation hence 
concretely depends on the parametrization chosen of the branch points in terms of the 
opening angles and the side lengths. Let us first consider the behaviour close to a branch 
point, say ujq. By a linear change of variables, uj = uj — ujq before the variation, and 
O = uj — uj'q after the variation, we can consider the branch point ujq to lie at the origin, 
and not to move under the variation. Taking a small corner around the inverse image of 
the branch point (in the 2-plane), and ignoring terms that vary slowly over the corner, 
we have 

25crcorner = log |^|^ = log - log |-^|^ = 

az auj duo 

5aQ log |ao2;| "0 + ^ 5aj log \uOj — uoq\^ — ^(1 — %)[ h c.c] (A2.3) 

In the bulk (i.e. away from the corners), we do not shift the cj-coordinates, and have 
instead: 

Y.Saj log \uj - u;,f - - aj)|^i^ + ccj (A2.4) 

■ UO — UOa 

3 3 J 

However, by Appendix 3, in piece-wise fiat domains, the only bulk contribution will 



arise from the normal derivative at the boundary of |A2.4| , so let us compute that (in the 
tu-plane, that is, at the real axis): 

2n . V5ab^ik = -2 E HlT^ + 2 E(l " ^^O^t?-;^] (^2-5) 



If uOj is on the real axis, [A2.5| vanishes identically. If uoj lies in the upper complex plane 



(an interior corner) we are to integrate the normal derivative along the boundary. The 
second sum in A2.5| then clearly gives zero, while the first integrates to 



/ 2n ■ V5(T]^^i]^(is = -47r ^ ^"j (A2.6) 
For simplices with disc topology, all the contributions except the rather simple result 



A2.6| therefore come from the corners, and for polygons, there are no interior corners, and 



[A2.3| gives the only relevant terms. We have not been able to put |A2.3| in a form which is 



manifestly a total variation, but in all the special cases where we checked it, it turns out 
to be so. In some cases it is possible to choose the parametrization of a family of polygons 
such that the corners on the boundary do not move: in this case all the contributions are 
from the first two terms in |A2.3| . For triangles we have the final simplification in that 



we can use the parametrization (|A2.2|) , and then the only contributions are from the first 
term in |A2.3| . 
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A3. Contribution of a straight boundary 



On piece-wise flat surfaces, away from the corners, we can always approximate tlie lieat 
kernel with the free space heat kernel, or the heat kernel close to a straight boundary. 
The error term {Kd — Ko in |A1.6| ) is exponentially small in |. In free space the diagonal 



elemeent of the heat kernel is Ko{x, x, t) = That leads to a contribution to Z[x, x, s) 
from Kf) in |A1.6| , which has a singularity at s = 1 with residue ^ , but which is regular at 



lower s. Inside the piece-wise flat surfaces there are therefore no contributions to Z(0), 
or to the variation of Z'{0). 

The heat kernel in the half plane with Dirichlet conditions on the x-axis is expressed 
as the difference between a "direct" and a "reflected" term: 

K{{xi,yi),{x2,y2),t) = ^(exp( )-exp( )) 



The diagonal elements are then 



(A3.1^ 



Kiix, y), (x, y), t) = ^(1 - exp(-^)) (A3.2) 

We compute the Mellin transform multiplied with a convergence factor exp(— /it), and s 
between ^ and 1: 



To investigate the contributions to Z{s) we may integrate y over a strip along the boundary 
with width 6 and length L. We then obtain: 

-L ^ -J - (A3.4) 

47rr(s) 2s- 1 ^ ^ 

When we continue to lower values of s we find a pole at s = | with residue — ^ which 
implies that the short term asymptotics of the heat kernel has one term which is a line 
density at the boundary, with coefficient —t~^-^. 

At the origin, the contribution goes as s, hence there is no contribution to Zd{0) from 
a straight edge. We obtain contributions to 2'^(0) as 

(A3.5) 



AttS 

The interpretation of A3. 5 is, that apart from the finite and regular terms in Z'{x,x,0), 



there is a distribution, which is defined by the analytic continuation of ?/^* ^. If we 
integrate over strips parallel to the boundary, the terms as in [A3.5| would cancel pairwise. 



and we would be left with the outer boundary term. Hence we conclude that a straight 
boundary gives a contribution to Z'^{x, x, 0) which is no more singular than a line density. 
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and hence gives a contribution proportional to boundary length. In general one could 
integrate Zd{{x, y), (x, y), s) against a smooth test-function /(x, y), with s in the interval 
[|, 1], and then analytically continue to smaller s. Such a test-function is the conformal 
variation 6(j{x,y). Assuming an expansion normal to the boundary at x as 26a{x,y) = 
doix) + y ■ ai{x) + . . . we find 

f dx f dy2Sa{x,y)Z{{x,y),{x,y),s) 
Jo Jo 

= ['dx f dyl-'^^^iyy-'jiaoix) + y ■ a,ix) + . . .] 
Jo Jo 47ri [s) 

Yd _ s) 52*^1 Yd - s) S'^" 

The contribution to the variation of Z'j^{Qi) is the limit of [A3 .61 as s tend to zero, which may 
clearly be expressed in the normal derivative at the boundary of the conformal variation: 

1 

'^^straight boundary(O) = ^ ' (A3.7) 

A4. The Sommerfeldt kernel 

A. Sommerfeldt in 1896 solved the problem of diffraction of light by a perfectly conducting 



half-plane [p!q] . The solution takes the form of a kernel periodic in the angle variable with 
periodicity 47r; the difference of one "direct" and one "reflected" wave vanishes at and 
27r. 

We will need the solution to the diffusion problem in a sector with opening angle 
vra, which is quite analogous, but for completeness given here. (The solution of diffusion 
problem at an interior corner with total angle 27ra is obtained in the same way, by keeping 
only the "direct" term.) If the opening angle is of the form - the sector can be refieced in 
its side 2n times to precisely cover 27r, and the solutions to both the diffusion problem and 
the diffraction problem are expressed by the method of images. Sommerfeldts solution is 
a substitute when the reflections do not make up a full turn, and is given by a certain 
finite number of image charges, and a correction term. It has the following integral 
representation: 

Ks{r, ct^; r\ 0'; t) = ^ exp(-^;ii^) - 0') - + 0')] (A4.1) 

where the important part is 

1 [ t jr. dS 

— / exp(acoso) 

27ra Ja+b 1 — e 



^"(^'0) = TTTT /. „exp(acos(5)- (A4.2) 



A and B are paths in the plane of complex 6 that go asymptotically from vr + i ■ 00 to 
—71 + i ■ 00, and — vr — z ■ 00 to vr — i ■ 00, respectively. Essentially this is a superposition of 
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free heat kernels between x and y', where \y'\ — r' and \x — y'\'^ = + r'^ — 2rr' cos (5). In 
the bands ^^^^ < \Re{6)\ < ^^^^ the contour integral can be taken to infinity since 
^{\x — y'p) oo, but not in between. 

Ks satisfies the heat equation because it is superposition of free heat kernels, it is 
symmetric in {x,y) and periodic in and (p' with period 27ra by construction of Ua, and 
it vanishes at the boundaries of the sector because it is the difference between a direct 
and a refiected term. Furthermore, away from the imaginary axis it is analytic in a. 

By deforming A and B into the straight lines tt + iy, — tt + iy and [— tt, tt], can be 
written as 



i/a(a, 0)= ^ exp(a cos(27rQ;A; — 0) 

k:-l<2ak-4>/Tr<l 

sin(7r/a) f°° cxp(— acosh?/) 



27rQ; J-oo cosh(?//tt — i(p/a) — cosir/a 



dy (A4.3) 



From the image charges of the "direct" term, it follows that the normalization of the 
kernel is correct. If for some k an image charge wanders through the line ±7r, we ought 
to take half the residue and the principal value of the integral in the usual fashion, but 
then a equals ^, and the prefactor of the integral is zero: that is; we have just a solution 
by images. 



A5. A useful integral 

In this appendix we discuss an elementary but useful integral. We would like to compute 
a contribution to .^d(O) or Z'j^iQ) by substituting the true (unknown) heat kernel in a 
compact domain, with essentially something like a heat kernel in the full plane. This 
includes the case with a flat boundary, where the heat kernel can be written as the 
difference between a free heat kernel and an image charge reflected over the boundary, 
and also the case in a sector, since the Sommerfeldt heat kernel can be written as a linear 
superposition of free heat kernels. The first problem is that the eigenvalues in an infinite 
domain accumulate to zero; the Mellin transforms are therefore not convergent for any 
parameter s. We are however interested in Mellin transforms of heat kernels in compact 
domains, which are asymptotic to free heat kernels for short times, but fall of exponentially 
for large times. Furthermore, we are only interested in the singular contribution, so the 
problem can be resolved by limiting the integral in the Mellin transform to the interval 
[0, 1]. In practice such an integral is not analytically tractable. Instead we can introduce 
a convergence factor e"''*, including which the MeUin transform converges for Re{s) > 1, 
and look for the singular part in the /^-independent part of the result. 
We therefore consider 

r f-^e-^-'^'dt = 2(-f^KJUaii)h = Us)/!-' + r(-s)(a)^ + /i) (A5.1) 

by the power series expansion of the modified Bessel function. 
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Two variants, with s around the origin, are 

i [s) Jo i [S) 

and 

i (s) JO t (■5 — 1) i (s) 

As fi tends to zero, in the first case, the value at zero is 1, and the derivative is 27 + log a 
by the expansion T{s) ~ ^ — 7 + . . . In the second case the value at zero vanishes, but the 
derivative is -. 



A6. The singular corner contribution 

We will here compute the contributions to Z£,{0) and 6Z'^{0) from a corner. For Z£,{0) 
we are by Appendix 2 and Appendix 4 looking for the limit as s tends to zero of 

/ cfxZ{x,x,s) (A6.1) 
Jcorner 

where the zeta function density is computed from the Sommerfeldt heat kernel, and the 
corner is delimited by a cut-off radius Tq. 



For SZ'jj{0) we use the variational formula (|4l|), and look for the term finite in s as s 
tends to zero of 

d^x[2Sa{x)Z{x,x,s)] (A6.2) 



i 



corner 



An inspection of |A2.3| shows that the conformal factor close to the corner contains two 
types of terms, one which is proportional to the logarithm of the distance to the tip of 
the corner, and one which is constant over the corner. 



We choose in the following to keep the constant factors in first term in |A2.3|, so that 



what we compute will directly be the stricly local contribution to 6Z'j^{0), which can then 
be given in integrated form. The contributions of the last two terms in |A2.3| will be 



proportional to the contribution to ^^^(O), and are easily included by comparison. 
Our computations therefore boil down to |A6.1| and 

i-TTa i-To ^(y 

/ / rdcpdrl — 2\ogaZa{x, x, s)] (A6.3) 
Jo Jo a 

rira pro 

/ / rd(j)dr\— log r^Z J x,x,s)] (A6.4) 
JO JO a 

where we can directly integrate over since our test functions (a constant and logr^) do 
not depend on 0. 

It will turn out that the reflected term in the Sommerfeldt kernel only contributes to 
6Z'j^{0), and only by a simple term. We can therefore in parallel compute the contributions 
to a corner at the boundary with opening angle vra, and a corner in the interior with 
opening angle 27ca, which almost only differ by a factor of two. By analyticity in a we 
can choose the convenient range ^ < a < 1. 
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The image charges 

The charge from the direct term, and the image charges from the reflected term give 
contributions to the MeUin transform, integrated over the area of the corner, as 

— / rdr dcpll- V exp(- — (1 - cos(27raA; - 20)))] (A6.5) 

The direct term charge can be neglected, as along a straight boundary. With | < « < 1 
the only relevant values of k in the reflected term are and 1. For k = the integral 
over (f) is limited between and |, while for A; = 1 it goes between ira — | and ira. With 
y — rsin(0) in the first case, and y — r sin(0 — ttq;) in the second, we find with s between 
I and 1: 



lJdcl>rZ^{{r,cl>),{r,cj>),s) = -^^^-^r'^-^ 2 J dcP sin''-' 



47rr(s) 

r(i - sy^-^rj^ns - 1) 

47rr(s)r(s) 



(A6.6) 



When we integrate up to a radius ro and continue to lower s we find, as we should, a pole 
at s = 1 with residue — 

To investigate the contributions close to zero, we use the formulae 



/ 

Jo 



/ 

Jo 



2s 2s 2 ^ ° ^ ' 
" drr'-Hogr' = |[|] = + i(logro^)^ + 0(s) 



where we understand that we first integrate over r, and then analytically continue to the 
neighbourhood of the origin. 

The contribution to 5Z'^{0) from the reflected image charges (hence only for boundary 
corners) is thus 



(A6.7) 



5a 

The reflected term: the integral 

The integral part of the reflected term comes from 

— ^ / rd(p / — , , , ^dy, (A6.8 

Sn'^at Jo J-oo cosh{y/a — 2i(p/a) — cosTi/a 

where we have integrated over the angle variable 0. The 0-integration forms part of 
closed contour from —ioo to the origin, along the real axis to ira, and down parallel to 
the imaginary axis to ira — ioo. The integrals over the two lines cancel. If y > there are 
two poles inside the contour, but their residues cancel. The angle integral is thus zero. 



r"^(l+coshj/) 
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The direct term: the integral 

We consider 



sinvr/a /™ exp(- - (1+-^.) ) 

' ' rd(j) — ^ ^ j-dy (A6.9) 



Sn'^at Jo J-oo coshy/a — cosn/a 

Following Appendix 2, the Mellin transform of [A6.9| , performed termwise inside the inte- 
gral over y, gives 



r 



2s- 



_iT{l — s) sin(7r/a) 



- s sm(7r/a) f°° i , / , ^ 

i (sjoTT J-oo (cosh y/a — COS TT/aj (cosh y/2j 



The integral is most tractable if considered as a correlation integral between the two 
factors in the denominator, then by Percival's formula^ we have 

2s-ir(l-s)a 2-2s P on , • 1 • w sinh Try cosh Tray -sinh Tray cosh Try 

1 (sjoTT J-oo smhavry 

(A6.11) 

This expression can be written more compactly as 

where the numerical values are determined by the integral 

, f°° ^, , ,sinh TT-u cosh 7r«v — sinh Tray cosh TT-u, , 

A{s) = B{l-s + ty,l-s- ty)[ ^ -]dy A6.13 

j-oo smh any 



Using the identity 



we have 



B{l + iy,l-zy) = -^ (A6.14) 
smh Try 



/oo 
[ny coth any — ny coth. Tcy]dy (A6.15) 
-oo 

and 

/oo 
(2^(2) — ^(1 + iy) — ^/'(l — iy))[ny coth. any — ny coth.ny]dy (A6.16) 
-oo 

where + 1) = -7 + 



^00 ^ i_ _ 1 ' 

^71 n+x' 



^This transformation brings us closer to the Lebedev-Kantorovich tranform of D. Ray, that was used 
by McKean and Singer in 
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Contribution to Zb{0) 

We integrate |A6.12| over the radial variable: 



(A6.17) 



A{0) can be computed by introducing a factor e**^^, and closing the integral over y in 
the upper half plane. The residues contribute 



OO -| 



n=l " 



TT 



2 a2sinh^(e/2a) sinh^(e/2)' 



which, in the limit as e — >^ 0, becomes — 1) 
The value of 2'q,(0) is therefore 
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— a) 



a result first derived by McKean and Singer |]T5 



(A6.18) 



(A6.19) 



Contribution to ^Zj^(O) 



Using again |A6.7| and the expansion of Ai^s) and the prefactors at the origin, we obtain: 



[— ][S7^^(0) + ^(^'(0) + A(0)(27 - 21og2))][--^ + ^] + h.o.t (A6.20) 



a '" 27r 27r 
We have a pole in s: 



2s2 



1 (5a a 

s a 47r 



A(0) 



(A6.21) 



Since ^A(O) is the contribution to .^_d(0) from the corner, we can write this variation as 



\ba , 1 ^ 



(A6.22) 



Let us work out that this variation is only seemingly in contradiction with the known 
expression of ZJ^^ in |A6.19| . The variation of 0) is of course a point mass at 

the corner, with weight |f (~~2 ^ !)• The variational forms and ( |5BD are however for 
the variation of Zd{^s\ not for the zeta function density. The variations we consider go 
between different simplices with disc-like topology, and therefore leave the sum X^i c ~ 
2aj)+Z]b c invariant (see equation ^). The analytical expressions of the variation 

are then only determined up to terms linear in the variations of the angles, which is the 
discrepancy we have in |A6.22. 
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The contribution finite in s in IA6.20 



IS 



Sa[^ (2 log a - 27 + 2 log 2) - ^] (A6.23) 

We compute A'{0) by introducing a convergence factor e*"^. Since ip{z) = ijj{z), it is 
enough to consider 

Cn = ( )\iTy cothany — Tiy coth7Ty]dy (A6.24) 

J-00 n n — iy 

in terms of which A'{0) = 2A{0) — 23?X]^i ^n- The two terms in Cn have residues in the 
upper half plane, and sum to 

I ^) 

2n^aHinh'^{e/2a) sinh2(e/2)^ 



00 



m( ^ g-W" I g-m^) 

a'^{n + m/a) {n + m) 



m=l 

= -f(A-l) + -(l-i) 

on a 

-27me^" e"'^" [— - , ^ J rf/x + 0(e) (A6.25) 

ie ^a(e'^/" - 1) (e^ - 1)^ ^ ' ^ ' 

With some care the sum in n is: 

00 

Y.Cn = -27r/(«) = 

n=l 

r°° 1 , 1 1 1,1 



^""io eA'-lf4sinh2(/x/2) 4^^ sinh2(^/2a) 12^^^ ^^^^^ ^^^''^^^ 
Now we can collect all the strictly local contributions to Z'jj{0) and write 

SZ'^iO) = + H^il^ - «)(log« - 7 + log2 - 1) - /(a)] (A6.27) 

Here we have taken a boundary corner, and hence included the contribution from the 
reflected image charges. 

It is clearly of interest to give [A6.27| in integrated form, whence we introduce the 
primitive function of I (a): 



00 



J(a)= I ^^[— coth(-^) + (A6.28) 

^ ' Jo e^-r2/x ^2a' 4sinh2(|) 12^a ^ ^ ' 

1(a) = — ;— ^(a) 
da 

We know from Appendix 3 that the contribution to Z'j^{0) {note: not the variation) of a 
straight boundary is proportional to the boundary length. The integrated form of [A6.27 
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should therefore be zero at a = 1. As in the variation of Zd{0), we should in addition 
expect undetermined linear terms in the angles. 

A rather long derivation, to be given in Appendix 8 as we investigate the special values 
of a = shows that we may integrate [A6.27| , compare with one of the integrable cases, 
include the integration constant in the contribution from the corners, and simplify to: 

Z'M = - «)(7 - log 2) - + 3 + a) hga + J{a) (A6.29) 

where the last term has the more symmetric integral representation 

roo \ \ jr X 1 1 

J(a)= / (coth(— )-acoth(-))--(--«)]rfx (A6.30) 

Jo — 1 2x 2a 2 12 a 

Difference at interior corners 

The only change is that the angle-independent direct term is integrated over twice as 
wide an angle, and the contribution from the reflected image charges are absent. The 
contribution to Zd{0) is then 

^a(0) = -^(--«) (A6.31) 
12 a 

and the strictly local contributions to Z'^{0) 

^a(O) = li- - «)(7 - log 2) - i(- + a) \oga + 2J(«) (A6.32) 
b a b a 

Summary of corner contributions 

For simplices with the disc topology, the contributions to the variation of the zeta function 
derivative are given by |A2.3| and a simple sum over the interior corners. There are three 
terms in |A2.3| , of which the first is the only one that arises in triangular domains. It is 



written in final variational form valid for corners on the boundary in |A6.27| , |A6.26| , and 



in integrated form in A6.29 , |A6.30| . The last two terms in |A2.3| are of the same type as 



A6.3| , i.e. in variation form they will give something proportional to the contribution to 



Z£){0) from the corner. We give those expressions here, under the (slight) simplifying 
assumption that none of the branch point in the cu-plane lies at infinity: 

•^^extra = E ^".(0)[5«, log - u^' - (1 - a,)[^-^^^^^ + c.c]] (A6.33) 

A7. Computation using the heat kernel 

This approach was used by DowkerQ] to compute the contribution to Z£,{0) from a corner. 
We begin with the asymptotic formulae for the trace of the heat kernel in a domain D: 

Tr4Knix,x,e)] ~ ^ + ^ + ^,,(0) + O(e^) (A7.1) 

e £2 
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Tr,[2Sa{x)KD{x,x,e)] ~ Ail + - log e 5[Zc(0)] - 75(^^(0)] 

e 62 

+ 5[Z^(0)]+O(6^) (A7.2) 

As in Appendix 6 we only need to consider the image charges from the reflected term and 
the integral from the direct term. 

The direct term: the integral 

It is convenient to go back to the original form of the Sommerfeldt kernel (including the 
hne integrals and the direct charge) : 

^direct(^' ^' ^) = / ^^e-w(i-^°^^) - (the pole at 5 = 0) (A7.3) 

We are to integrate the kernel, by itself or multiphed by ^ loga^r^, over the corner, and 
then catch the finite piece as t tends to zero. The integration contour is chosen such that 
the prefactor of in the exponent always has negative real part. We can therefore extend 
the integral over all r (as in Appendix 5), without introducing significant errors. For the 
contribution to Zd{0) we then find 



7ra 



2 



racp I K^:^ Jx,x,t) = -— / + Oie-"^- * (A7.4) 

Jo airectv / Svr Ja'+b' l _ « 1 - cos 5 ^ ^ ^ 

where we have afterwards moved over the integration contour to the straight lines parallel 
to the imaginary axis at tt and — tt. Except for an error exponentially small in |, we 
have only left a term independent of t, that will give us what we want. By a change 
of variables, S — > —5, the integration contour A' goes into —B', and vice versa, and the 
second factor in the integrand ( i_cosg ) is unchanged. It is therefore sufficient to keep in 
the first factor, the part odd under this reflection, and we have 

-1 A 1 

dS cot{—)- (A7.5) 



Idni Ja'+b' 2a 1 — cos 5 

This integral is now in fact convergent at infinity, so we can close it, and evaluate it by 
computing the residue of the pole at the origin. Hence 

pira pro 1 S 1 

/ rd(f) / drKj,: Jx,x,t) = -— 27riRes[cot(— )- -]s=o 

Jo Jo uiit;L.iv 2a 1 — cosd) 

which is, of course, the right answer. 

To compute the variation of Z'jj{0), it is again convenient to extend the integral over 
the radius to infinity. We then find 

rd<Pr r -log{a'r')e-'^^'-^"^'^ = , , ("7 + log , ^""'^ J (A7.7) 

Jo Jo a 1 — cos 1 — cos d 
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(we do not write out the variational factor 5a, except where necessary), and hence 



7ra rr^^ 



rc 







dr[26a{x)K^:^^^^^{x, x, t)] 



/ ^ + A7.8 

Sira JA'+B' 1 — 1 ~ coso 1 — coso 

PuUing out the term divergent as logt we have 

5a 1 , 1 , / . » „N 

logt— - - -a A7.9 
a 24 a 

which is equal to — \ogt 6[Za{0)], up to the term linear in the the angle. 



Taking into account that the finite part of |A7.8| is 5[Z^(0)] — 75[Za(0)] we have 



5[Z;(0)] = 5[Z„(0)](-log(2a^)+27) 

6a f d6 1 , 1 , . „ X 

log- (A7.10) 



87ra JA'+B' 1 _ g-^r 1 — cos 5 1 — cos 5 

We choose the branch of log symmetric under 5 —5. Then we can again substitute 

for — ^"Tj its anti-symmetric part, and the integral is convergent at infinity. Essentially 

1— e a 

log i_losS ^ function of 5^, so it can be chosen to have a branch cut along the negative 
real axis in the 5^-plane. That translates to two branch cuts in the 5-plane, along the 
positive and negative real axis. Evenness under 5 — > —5 then determines the phase choice 
in the left half plane. The best we can do now is to pull both integration contours in to 
the imaginary axis and integrate over the branch cuts. In addition we have a singularity 
at the origin. We separate out a small circle around the origin with radius e, and an 
integral along the branch cuts from ±ie to ±200. The integrals along the branch cuts give 

1 re- / ^ ^ ^ 

— I do cot — 



4a Jit 2a 1 — cos 6 

Close to the origin we estimate {6 = ee*^): 

log — - = Iog2-21oge-2z0 + — + ... m>0 

1 — cosd 12 

= log2-21oge-2z^ + 2z7r + -3-5^ + ... m<0 



12 



which eventually gives 



— (^-l)(log2-21oge)- — + At (A7.12) 

24 a2 ^ ^ ^ 24 2e2 ^ ' 



28 



Collecting all terms from the direct integral we have 

•^[^alO)] direct integral = "^"[^(^ " ^)(~27 + 21og(2a) - 2 log e) 

2e2 ^ 4a 12 ^a^ ^ 



+Ar''.-th(|)-^| (A7^13) 



5a^ 



The image charges 

We should compute the integral 

'^[^:.(0)]image charges = -^"t// ^d<p rfr— log(ar)e-ir(i— 2,)] (^7.14) 

The problem is that we have to separate out a rg-dependent piece diverging as t"^, before 
we can have the finite value at the origin, and as it stands it is not possible to extend the 
integral over r to infinity. If the integral would have been over a small square, we could 
have separated in x and and then the integral would be much easier. Let us first see 
that the segment of the square outside the quadrant is unimportant, and then compute 
the integral over the square. The difference is 

-/ dy dx—-log{a\x^ + y^))e-'r (A7.15) 

Jo J^/rl-y^ znat 

which is essentially 

- r dy-^[\og{a^rl) + ^ + . . .je"^ (A7.16) 
Jo Anatro 

which only gives a contribution of order t^. 
The integral over the square is 

-/ dx dy- log(a^(x^ + ?/^)e"^ (A7.17) 

Jo Jo znat 

After first integrating over x, we can extend the integral over y to infinity and find: 

- / dx (loga^ + 2ro(logro - 1) + ny + . . .)e~~ = Const. -T^ (A7.18) 

JO 27rat 4a 

Summary 

The terms diverging with the cutoff e can be brought into the integral, and then give, 
after some algebra, 

6[ZM] = - 1)(1 + 7 - log2a) - ^ + J'{a)] (A7.19) 



which is the same result as in |A6.27 
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A8. Corners with opening angle ^ 

We proceed to simplify the the general expressions of the contributions to the regularized 
determinants, when the angles are of the form -. 

We recall that the contribution to Z'j^{0) from a (boundary) corner is 

11 1 

A + Ba — - log a H (7 — log 2 — log a) a log a + J {a) 

4 12a 12 

where the last term has the integral representation 

J(a)= r — ^ — [— coth(— ) - —i- + a)\dLi (A8.1) 

^ ' Jo e^'-V2fi ^2a' 4sinh2(|) 12^a ^ ' 

and have included an undetermined term linear in the angles {Ba), and an integration 
constant, that we will not need explicitly {A). 

The analysis proceeds by computing the integral for a being -. Let y stand for e'^. 
Decomposing various rational functions of y: 

n-l 

log(?/" - 1) = 5Z log(l/ - A^) = 

n ^ 1 
r - 1 h A.I/ - 1 
I , n n, 1 1 "1^1 A„ 1 

+ ^) = -, TTTT + ^r^ TT-Er 



y-X^yn_X 2' 2(y-l) ^^\-\y-\ 



A, ^ 1-n 
.=il-A. " 2 

we can absorb the terms proportional to ^ into 5, and write the modified integral: 



n n 



- — y — — ^ — 7 — ^]dx A8.2 

nJe ^ X Ky- \, 12 y-iy ^ ^ 



The cut-off e is introduced for later convenience. In fact, we can go backwards and express 
imias 



Since up to simple terms in a that we keep, |A8.3| and [A8.1| agree on the integers, and both 
expressions are analytic in ^, then they must agree overall. We are therefore allowed to 
substitute for |A8.1| the more symmetric expression |A8.3| . Continuing the analysis of |A8.2 
we use 

r'^ dx Xy f'^ dx ^ ^ e^^"")^ 

Je X y — Xy Jne X ~^ 6^ — 1 



30 



oo g(l-a)a; 1,1 1 

'- — a)(7 + loge) + logr(a) — - log(27r) + small as e 



le X -1 e '2 ° ' ° ^ ^ 2 

A™ , 1+n 



and express the integral J as 

v=l p 

which may finally be simplified to 



Absorbing further terms proportional to ^ into S, we have finally 

Zi (0) = A + 5- + -— log 2 + (- - — ) logn + E (- - ^) log r(^) (A8.6) 
n 12 4 12r2 ^ 2 n n 

The integration constant 

We need one exact value to fix an integration constant. For the equlateral triangle, the 
regularized determinant is (see Appendix 11): 

^Equilateral (0) = ^ log - ^ log 2 + ^ log 3 + ^ log ^ (A8.7) 



If we on the other hand use |A8.6|, we find 



^;quilateral(0) = H (0) = 2 + ^ log 3 + 1 log ^ + 3i + 5 (A8.8) 



^3 

SO we have the integration constant 

3i + 5 = ilog2 + ilog7r (A8.9) 
3 2 

Simplified general expression 

We choose to include the integration constant in the stricly local contribution to Z^(0), 
by adding 

(l-a)[^log2 + ^log7r] (A8.10) 
o 4 
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to a boundary corner (and twice this quantity to an interior corner). The final expression 
with an opening angle ^ then reads 

11 1 1 — n 11 

= ^^' n^^e + 4 + + ^4 - 

+ E(^--)iogr(-) (A8.li) 

^ 2 n n 



For a general opening angle we go back to |A8.3| , keep track of the various terms propor- 



tional to - that had been absorbed into B, and arrive at 

KiO) = - «)(7 - log2) - + 3 + «) log« + J» (A8.12) 

12 a 12 a 

Throughout this appendix, we have simplified on the integral J, which is the same for 
corners on the boundary and in the interior. The strictly local contribution to Z'^{0) from 
a corner in the interior is thus 

= li- - a)(7 - log 2) - li- + a) \oga + 2J{a) (A8.13) 

A9. Asymptotics of ^^(0) 

In this section we investigate Z'^{0) when a is large or small or close to one. 
Large and small a 

It is convenient to divide up the terms in Z'^{0) according to whether they are symmetric 
or antisymmetric under the transformation a — -. If we take the integral J, it makes 
sense to first introduce a finite symmetric alternative expression 

dy 1 1 1 , I ^,1 1 , 1 , , , . ^ 

"^^^ = / —^JJ^— — + i^ + V^)- + -ia + 3 + -)\oge{A9.1) 

Je y ey^^ — 1 _ X 2e^ y/a 2e 12 a 

where the limit as e ^ is understood. Using the elementary integrals 

dx 



- 1 



loge + 0(e) 



dx I 11, 7 1 , x 

+ -loge+ ^ - -log27r + C(e) 



e^-lx e2 2 2 

oo 



dx 1 1 1 1 , 1 7 ^// X N 

- - loge + +- - - C'(-l) + 0{e) 



e^-lx2 2e2 2e 12 ^ 12 12 
dx 11 1 7 X X 

(A9.2) 
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we can express the difference between J and J5 as 

1 1 1 + a ^ 

J{a) = Js{a) + (a + 3 + -) loga + ^(7 - log27r) - + C'(-l)) (A9.3) 

and divide up Z'^{0) in symmetric, antisymmetric, and linear parts: 

Z'^iO) = [J^(a) + ^(- + 3 + «)-^log27r] 
VI a 4 

- «) log2 + + 3 + a) log a] 

-a[ilog27r + C'(-l)] (A9.4) 



The behaviour of J5 as a turns to zero is found expanding the integral in [A9.1| , and 
subtracting the terms in the expansion divergent with e: 

Mcy) -a^o ^(^ + 3 + «)iog« + (^-C'(-i))^ 

+ i(log27r-7) + E^7^«" (A9.5) 
4 ^3 n{n + 1) 

and from this we find 

ZUO) ~a-.o -d(l-log2)-a-l)) + a(^^±l^-^log27r-C'(-l)) 
a VZ 12 4 

V ^M^a" (A9.6) 
Using the symmetry under a — > ^ we have the asymptotic expansion for large a: 

Z'M ~™ --^(- + 3 + «)log« + «(i-(l-log2)-^log27r-2C'(-l)) 
VZ a VZ 4 

« 12 ^3 n(n + 1) 

Development of Z[^^{0) 

It is convenient to go over to the first definition of the integral J{a) and write 

KiO) = -^(- - a)(7 - log 2) - + 3 + a) loga + J(a) - aAJ 

VZ a VZ a 

where J has the integral representation 

T/ \ r°° I r I w X « 1,1 , 

J(a) = / — cothi — ) 7T—— 1 — \- a)\aa 

^ ^ Jo e^'-V2fi ^2a' 4sinh2(f) 12 ^ 
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and the difference has the integral representation 
AJ = / dx[-{— — — + 



x'(e^ -l)2 2(e^ -l)' (e^--l)3 Q{e^ - 1)' 
= -i7-^ + ^log(2vr) + C'(-l) (A9.8) 

The derivative of J(a) can be expanded around a = 1, and the successive terms evaluated 
in Mathematica, which gives 

J'(l + 6) = -l6 + le2 + 0(e3) (A9.9) 

and putting the various terms together we have for a corner on the boundary: 

ZUM = (^log2-^-llog(2vr)-C'(-l))6 

+ (- + + (-^ - Tz}2^y + 0{e^) (A9.10) 

^72 12 ^ ^ 144 12 ^ ^ ^ ^ ' 

The expansion for a corner in the interior is easily obtained by adding the expansion of 
I log(l + e), and doubling that result. 



AlO. Special cases 
Triangles 

For triangles one can choose a convenient parametrization in terms of the Schwarz- 
Christoffel transformation that maps a point {z) in the triangle to a point (uj) in the 
upper complex plane: 

^ = cji-"«(l -cj)i-°i (AlO.l) 
This fixes the area in terms of the angles to be 

Area{a,, a^, a,) = -— Vrn T ^^^'^ 

2 r(l - ai)r(l - a2)r(l - as) 

One corner (with opening angle Trao) will map to the origin in the cj-plane, one corner 
(vrai) to one, and the last corner maps to infinity. The parametrization is uniform over 
the space of triangles, in the sense that the branchpoints do not move. Hence the third 
term in |A2.3| is identically zero. The logarithm of the distance between the origin and one 



is zero, hence the second terms in [A2.3| are also zero. All we have left is then the strictly 



local contribution, so for this choice of the normalized areas, we have for triangles: 

p=l,2,3 



34 



An integration constant was included in the definition of Z'^{0) in Appendix 8. 

By numerically solving for about the first one thousand eigenvalues of the Laplacian 
in isosceles triangles, and then estimating the analytic continuation of the zeta functions, 
Luck found that the quotient 

Ct = ^ (A10.4) 



varies surprisingly little over the triangles (still taking the area A-1(J.2| ) P!1[ . We can express 
this Ct as 

Ct = ^"^''^''5''°' (A10.5) 

It has a maximum for the equilateral triangle, where it equals 4.591151 . . . 

The minimum of C,t is obtained as an angle tends to zero, and the value follows from 
the asymptotic expansion [A9.6| : 

Limit„_o^^ = 2(1 - log 2) - 24C'(-1) = 4.583813 . . . (AlO.6) 

Polygons 

In this case it is convenient to take a parametrization where the interior of the polygon 
[z] is mapped onto the interor of a circle [u] by a transformation that satisfies 

^ = TT(M-e''^'')i-"" (A10.7) 
dz V 

where e**^" are the branchpoints, and the interior angles are Tra^,. The important special 
case is a regular n-polygon, P„, that can be mapped to the unit circle by a transformation 
that satisfies 

_ = rr M-e— - = K - 1 " A10.8 
This parametrization fixes the radius of the circumscribed circle: 

r(i + i)r(i-^) 

Rn = ^ (A10.9) 

and the area, which is proportional to R^. We can now choose a family of polygons that 
smoothly interpolate between a regular m-polygon (at a = 0), and a regular mn-polygon 
(at a = 1): 

^ = (m'- - 1)1(1-'^) (m"^" - 1)^" (AlO.lO) 

which is convenient, since the parametrization is then uniform in the family. The variation 
of 2'^(0) will therefore be determined by the strictly local terms, and by the second term 
in [A2.3| . Let us express 2'^(0) as a function of the parameter a in the family, and write out 
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d rrl 

iWllig, UUiiLllUULlUli LU 

from the comer which is mapped to the branchpoint e 



these two terms in |A2.3| in more detaih they give the following contribution to ^-^z)(a)(0) 



lT;iv 
mn ' 



^^z)(a) (O)]from corner v " ^^"-('^)*^'^'' ^ ^ 



log M„/ - 

uu , , da 



which can be rewritten 



^^D(a)(0)]from corner v " ^[^"4")*^^) + 12*-^^^ ~ a„(a))A^(a)] 

6 da ^ 
We have to differ between whether n does or does not divide v. For the first case 



Xy{a) = - — (l-a)log 7-\^=i a log 



2 


:i- 


a) log 


m 






2 




a 

a H — 


m 




n 




— a] 


log e 



|m=1 



u — 1 mn u — 1 

2a , , / . X 
logn n f (A10.12) 

mn 

while for the second 

A^(a) = — — (1 — a) log |e^ — 11 ^logmn n Vv (AlO.13) 

m mn 



The expression [A10.11| contains one part which is a total variation, and an extra term. 
Let us first do the second, summed over the corners v. 

A^(a)^^^^^^ = ^^(- — (l-a)logm-— logmn)(— (1--) 
b ^ da b , m mn m n 

^ n\v 

+ ^ V(- — (l-a)log|e^ - 1| - — logmn)( — —) 
b ^ m mn mn 

which can eventually be simplified to 

21-a, 21-a, 

log m H log mn 

3 m 3 mn 

from which we have 

/' da]: E A.(a)^^ = hJ-iogmn-- logm) (AlO.14) 
Jo b ^ da 6 mn m 

The total variation in [A10.11| is summed over the corners, and gives in integrated form: 

2(n — I) 

nZ'JO) \ V log^ + Const. (AlO.15) 

3n{n — 2) 
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We are now ready to write down Zp^{0) for a regular polygon with n corners and radius 
of circumscribed circle R: 



Rl 3(n - 2; 



(A10.16) 



where Zp^iQi) 



n-l 



6(n-2) ^"^^ normal radius is as in |A10.9| . The integration constant 
turns out to be zero, since it has already been incorporated in the definition of Z'^iQ). 



Square and disc 

For the square we obtain 



- log(^l^) - - log 4 + 4(- log TT + — log 2) 
2 v^r(i)^ 6 ^ ^8 ^ 24 ^ ^ 

1 r(-) 1 15 

2 log ^ly + 4 log Area + - log TT + - log 2 



which agrees with the exact result from [A11.3 . 

The disc is obtained as the limit when n tends to infinity. We then find: 

z^jo) = i,ogis-2Az;(o)|„., 

= ilogfl + -| + ilog,r + ilog2 + 2C'(-l) 

6 12 I 



(A10.17) 



(A10.18) 



which agrees with [All.lO 



All. Integrable domains. 



In this appendix we collect the cases known to us, where one can directly deduce the 
derivative of the zeta function at zero. 

In a rectangle with side lengths A and B, the eigenvalues of the Laplacian with Dirich- 
let boundary conditions are 



2 9 
,m m 



^mn = 7r'{— + —) (All.l) 



The zeta function around the origin is 
^rectangle 

(0)= \ 
4ectangle(0) = i ' log[2-(f 

where rj is the modular form of Dedekind: 

00 

ry(g) = g2l J](l _ g"^) q = e~^^V f (A11.2) 

m=l 
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For the square, we have the simpler expression 

4quare(0) = \\ogA' + \ \og[^2^^^] (All.3) 

Three triangles tile the plane by reflections in the side: the equilateral (|, |, |); the 
bisected equilateral |, |, |, and the right angle isosceles (|, |, |)- In these domains one 
can solve for the eigenmodes of the Laplacian by superposition of plane waves ^ , We 
normalize the areas in terms of the side lengths (a) of the equilateral, the sidelengths of 
the legs in the right angle isosceles, and the sidelength of the longest side in the bisected 
equilateral: 

^(i,i,i),n,m = {—f{'n? + m'^-nm) n>m>0 
^(1,1,1),,,^ = {^f{n' + m') n>m>0 
^(1,1,1),^,^ = {—f{n^ + w? + nm) n>m>0 

The corresponding zeta functions can be written in terms of the Riemann zeta-function 
and Dirichlet L-series|TO[], which may in turn be resolved into sums of zeta functions of 
Hurwitz H{x, s) with different arguments x. Using 

OO 1 

H{x,s) = Y. 



H{x, 0) = l-x ^H{x, s)l=o = log 

2 as v27r 



and the normal areas determined by the representation of the Schwarz-Christoffel trans- 
formation; 

one finds after some algebra the following expressions for the determinants: 

%..,(0) = llog^ + log£%M?Z (AU,5) 

^3'3'3> S A(l 1 In Fa f -) 

%..,(0) = i,og-j^ + log£!iS^ (A11.6) 

5 , Area , r(i)3 24 2 9 7r2 

= 12'°«a7TT: + '°« 'r(i) (^"■^) 

I 2 ' 3 ' 6 ^ 



The eigenvalues of the Laplacian in on the upper half sphere of radius one are + If 
one imposes Dirichlet boundary conditions at the equator, one selects out the eigenmodes 
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odd under reflection in a plane through the equator. Each eigenvalue will then be / times 
degenerate. The zeta function of this spectrum is 



■^hemisphere 



(A11.8) 



from which one can derive pO| 



^hemisphere(0) = 2C'(-l) + 2log2vr-4 



(A11.9) 



The hemisphere can be mapped conformally to a disc with radius R. The difference of 
the regularized determinants on the hemisphere and on the disc can then be evaluated by 
computing the Liouville action (|l^) of the conformal factor [pO(| : 



^disc(O) 



^ log 2 + 1 log vr + 1 log + 2C'(-1) + 
bio 12 



(All. 10) 



Equation [A11.1(J| has also been checked numerically to great accuracy by computing the 
eigenvalues from the zeros of Bessel functions, and directly investigating the analytical 
continuation of the zeta function pi|. 
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